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Excitation of g modes in Wolf-Rayet stars by a deep 
opacity bump 
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ABSTRACT 

We examine the stability of I = 1 and i = 2 g modes in a pair of nitrogen-rich Wolf- 
Rayet stellar models characterized by differing hydrogen abundances. We find that 
modes with intermediate radial orders are destabilized by a k mechanism operating on 
an opacity bump at an envelope temperature logT w 6.25. This 'deep opacity bump' is 
due primarily to L-shell bound-free transitions of iron. Periods of the unstable modes 
span ~ 11 — 21 hr in the model containing some hydrogen, and ~ 3 — 12 hr in the 
hydrogen-depleted model. Based on the latter finding, we suggest that self-excited g 
m odes may be the sour ce of the 9.8 hr-periodic variation of WR 123 recently reported 
bv lLefevre etaD «200l . 

Key words: stars: oscillations - stars: Wolf-Rayet - stars: winds - stars: variables: 
other - stars: individual: HD 177230 (WR 123) 



1 INTRODUCTION 

Of all of the differing excitation mechanisms responsible for 
free oscillations of stars (for a comprehensive review, see 
lUnno et al.ll989l . and references therein) , none is more ubiq- 
uitous than the opacity, or 'ft', mechanism. Put in simple 
terms, this mechanism works in layers of a star's envelope 
where the opacity k increases in response to compression 
of the stellar material. Such an increase (initially produced 
for instance by stochastic perturbations to the equilibrium 
state) retards the escape of energy from the stellar core, 
wi th the trapped flu x being deposited in the layer as heat. 
As lEddingtonl ((1926,) discusses, the supply of additional heat 
during compression creates a Carnot-cycle heat engine capa- 
ble of converting part of the outflowing radiant energy into 
mechanical energy. In tandem with a restoring force that re- 
turns the layer back toward its equilibrium state (pressure 
for p modes, buoyancy for g modes) , the heat engine creates 
a condition of pulsational instability (overstability) . 

The role played by the k mechanis m in the classica l (S) 
Cephei pulsators was first identified by Zhevakin (1953). In 
these stars, the source of the instability is a peak in the 
opacity and its partial derivative kt = (9hiK/91nT) p at 
an envelope temperature logT ~ 4.65 where helium under- 
goes its second ionization. At lower luminosities, the same 
mechanism is primarily responsible f or the RR Lyrae and 
S Scuti classes of pulsating star (see ICoxlll98Ct and refer- 
ences therein). However, in the early- type f3 Cephei stars 
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and slowly pulsating (SPB) stars, the helium opacity bump 
is situated too close to the stellar surface to play any signif- 
icant role in destabilizing pulsation modes. 

In fact, the excitation mechanism responsible for these 
two classes of star remained a mystery until the ad- 
vent of new opac i ties f r om the OPAL and OP projects 
jRoeers fc Iglesiasl Il992t ISeaton et all Il994h . The treat- 
ment of spin-orbit splitting in same M-shell transitions 
of iron and nickel, absent in previous calculations, intro- 
du ced a new opacity peak at logT ~ 5.3. Stability analyse s 
by ICox et all dl992t). iDziembowski fc Pamvatnvkhl il993f) . 
and ^Dz!enibowski^taiT1^93Ti ' revealed that this so-called 
'iron bump' can destabilize p modes in j3 Cephei stars 
and g modes in SPB stars. More-recent investig ations by 
ICharpinet et all <1997l) and iFont aine et all (120031) have in- 
ferred that the same iron bump can also drive the pulsation 
of the short-period (EC14026) and long-period subdwarf B 
(sdB) stars. 

In parallel with these developments, the new opacity 
data also impacted understanding of the pre-white dwarf 
GW Vir (pulsating PG1159) stars. IStarrfield et"ai1 dl98l 
1984) found that g modes in these objects could be k- 
mechanism excited by an opacity bump around logT ^6.2, 
arising from K-shell photoionization of carbon and oxy- 
gen; however, an almost-pur e C/O mixture was re quired 
for the instability to operate ( Stanehcllini et al. 1991). This 
abundance restriction is relaxed when upda ted OPAL/OP 
date are adopted (see iGautschv et al"l l2005. and references 
therein), due to the appearance of extra opacity around 
logT w 6.3. The origin of this additional opacity appears 
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to b e spin-orbit effects in L- shell bound-free transitions of 
iron jRoeers fc Iglesiadll992h . 

From this brief survey of the literature, a natural ques- 
tion emerges. The iron bump in the OPAL/OP opacities 
leads to instability strips at both high masses (/3 Cephei and 
SPB stars) and low masses (sdB stars). Given that opacity 
bump at logT ~ 6 — 6.3 drives pulsation in the low-mass 
GW Vir stars, is there by analogy a high-mass group of 
stars that also exhibit pulsation driven by this 'deep opacity 
bump' (DOB)? In order for the pulsation to be remain non- 
adiabatic around the driving zone, and therefore allow the k- 
mechanism heat engine to operate, these putative stars must 
be exceedingly hot, with surface temperatures approaching 
10 5 K. In fact, these stars can be recognized as massive, pop- 
ulation I, core helium burning Wolf-Rayet stars. 

In this letter, we conduct a preliminary stability anal- 
ysis of a pair of model Wolf-Rayet (W-R) stars, to evaluate 
the hypothesis that g modes can be excited in these stars by 
the DOB. We introduce the models in Sec. |2] and describe 
the stability analysis in Sec. [3] Our results are presented in 
Sec. 2] and then discussed and summarized in Sec. |K| 



2 STELLAR MODELS 

We evolve a star with an initial mass M = 100 Mq and 
metalicity Z — 0.02 using a stellar structur e code written 
by JM . The code is described in detail in Ijimenez et alJ 
(2004); the details pertinent to the present work are that 
(i) opacities are interpolat ed in the revised OPAL tables 
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(ii) the heavy-element mixture 
is assumed for the initial com- 
position, and (iii) mass loss due to radiative driving is ac- 
counted for v ia a modified form of the formula given by 
lAbbotd dl982l) . The type-2 OPAL tables used by the code 
allow variation of the abundances of carbon and oxygen (in 
addition to hydrogen and helium), to reflect the effects of 
nucleosynthesis. Opacity changes due to a varying nitro- 
gen abundance are simulated by adjusting the carbon and 
oxygen values passed to the opacity interpolation routine. 
Since the CNO elements are minor (and approximately in- 
terchangeable) opacity sources, this procedure should not 
introduce any significant error. 

We track the star's evolution off the main sequence, 
across to the red region of the Hertzsprung-Russell (HR) 
diagram, and back again to the blue. Two particular evolu- 
tionary stages are then selected as exemplars of nitrogen- 
enriched (WN) Wolf-Rayet stars. The "WNL" (late-type) 
model is somewhat hydrogen-depleted, with a surface mass 
fraction X s = 0.117. The more-evolved "WNE" (early-type) 
model is completely devoid of hydrogen, but has yet to show 
any photospheric signature of carbon enrichment. Here, we 
use quotation marks around our models' labels to em pha- 
size that, in the same fashion as iLanger et all dl994l) . the 
designations 'early' and 'late' are to be understood in an 
evolutionary rather than a spectroscopic sense. 

The fundamental parameters of both stellar models are 
summarized in Table In Fig. we plot the opacity and 
its derivative kt as a function of temperature throughout 
the model envelopes. The DOB can be seen clearly as the 
peak around log T ~ 6.25, while the iron bump also ap- 
pears as a pronounced peak in the superficial layers around 



Table 1. Fundamental parameters of the stellar models intro- 
duced in Sec. [5] Note that the radii R* are for the hydrostatic 
core, and do not include t he extended wind re gion containing the 
t = 2/3 photosphere (see lBaschek et alJll99jl . 

model age/Myr logL»/L Af*/M Q R*/R Q X B 



"WNL" 
"WNE" 



3.57 
3.74 



5.77 
5.62 



20.6 
16.2 



10.3 
2.10 



0.117 
0.000 
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Figure 1. The opacity ft (top) and its temperature derivative 
= {9 In k/8 In T) p (bottom), plotted as a function of envelope 
temperature for the "WNE" (solid) and "WNL" (dotted) stellar 
models. The labels indicate the location of the deep opacity bump 
(DOB) and iron opacity bump (Fe). 



logT ~ 5.15. By exploring the effects of modifying the 
heavy-element mixture, we have determined that the DOB 
is due primarily to iron, with minor contributions coming 
from nickel and the CNO elements. 



3 STABILITY ANALYSIS 

To analyze the stability of the W-R stellar models in- 
troduced above, we employ the BOOJUM non-radial, non- 
adiabatic pulsation code developed by RHDT. The current 
revision incorporat es a few improvem ents over the baseline 
version described in lTownsen dl (l2005allrl . Most significantly, 
the roots of the characteristic equation 



V{uj) = 0, 



(1) 



that define the stellar eigenfrequencies, are now first iso- 
lated using an app roac h similar to those sugge sted by 
iDziemblowskil dl977l) and lShibahashi fc Osakl Jl98 lib In this 
equation, uj is the pulsation frequency normalized by the in- 
verse dynamical timescale rZL = (GM*/iZj) 1//2 , and T>(u>) is 
the dis criminant function defined by eqn. (18) of lTownsendl 
ll2005al) . 
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Figure 2. Normalized growth rates rj for the i = 1 (circles) and 1 = 2 (triangles) modes of the "WNL" (left) and "WNL" (right) stellar 
models, plotted as a function of pulsation period. The dotted line indicates the division between stable modes (rj < 0; open symbols) 
and unstable modes (rj > 0; filled symbols). 



The root isolation proceeds by dividing the complex-oj 
plane into many small rectangles, and evaluating for each 
one the contour integral 



C 



2TTiJ 



V(u>) 



(2) 



counterclockwise around the perimeter, via a recursively re- 
fined trapezoidal quadrature scheme. By Cauchy's theorem, 
C converges to the number of roots of the characteristic equa- 
tion enclosed by the rectangle. Once a root is thus isolated, 
the frequency 



1 



OJC = 



2m C 



lu AV(uj) 



(3) 



is used as the starting point in the secant root-finding al- 
gorithm implemented in BOOJUM. The C term in the de- 
nominator of this expression ensures that u>q lies within the 
enclosing rectangle. Even though this term has an analyti- 
cal value of unity, for a single enclosed root, its numerical 
value differs from unity due to the finite accuracy of our 
quadrature scheme. 

This root isolation procedure is extremely useful in stars 
having high luminosity-to-mass ratios (such as W-R and 
GW Vir star s), because the p ulsation tends to be so non- 
adiabatic fsee llJnno et alll989l their §22.1) that solutions to 
the adiabatic pulsation equations furnish very poor starting 
points for the secant root finder. However, the contour in- 
tegration often requires large numbers of quadrature points 
before C converges satisfactorily toward an integer value. 
Hence, although adding robustness to the process of solving 
the non-adiabatic pulsation equations, root isolation is prac- 
tical only when - as in the present work - a few individual 
stellar models are under study. 

For the "WNE" and "WNL" models, we focus the sta- 
bility analysis on g modes having harmonic degrees 1=1 



and 1 = 2, searching for modes with radial order s h between, 
appro ximately, and -25. Here, n is defined bv lUnno et alJ 
il989i t heir eqn. 17.5), within the gener alization to the 
ICowlin j jl94l|) nomenclature introduced bv lScuflairel (ll9T4t> 
and|Osaki ( 1975). Negative values of h typically correspond 
to g modes, but - as is often the case with more-evolved stel- 
lar models - n does not vary monotonically from one mode 
to the next, nor is it a unique index. 



4 RESULTS 

The results from the stability analysis are presented in 
Fig. [5] which plots the normalized growth rate rj = 
— $J(a>)/3i(a>) of the modes found against their period P. Ev- 
idently, the "WNE" model is unstable (rj > 0) against 1 = 1 
and I = 2 g modes having periods in the range ~ 4 — 12 hr 
and ~ 3 — 8 hr, respectively. Similar instability is seen in the 
"WNL" model for 1=1 modes spanning the period range 
~ 11 — 21 hr; however, the 1 = 2 modes of this model are all 
found to be stable. 

To explore the origin of the instability, Fig. plots the 
differential work for the 1=1 modes of each model that 
exhibi t the largest growth rates. As discussed bv lUnno et alJ 
(1989, their §26.2), the k mechanism is operative wherever 
dnr/dr > 0. From Fig. we see that this condition is met 
in the vicinity of the deep opacity bump, at a temperature 
logT « 6.3. Since this location coincides with the sharp 
positive peak in both models' differential work (indicating 
strong driving), we conclude that the unstable 1 = 1 and 
1 = 2 modes are K-mechanism excited by the DOB. 

In Table|5J we summarize the properties of the unstable 
modes. The long-period limit of the instability arises because 
the pulsation interior to the opacity bump, at a temperature 
logT « 6.5, becomes sufficiently non-adiabatic to activate 
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strong radiative damping there. Likewise, the short-period 
limit comes about because low-order modes exhibit a rela- 
tive Lagrangian pressure perturbation 8p/p that is too small 
in the vi cinity of the DOB to gene rate appreciable excita- 
tion (see iDziembowski et alJ Il993l . for a discussion of the 
role played by Sp/p). The £ = 2 modes of the "WNL" model 
are all stable for much the same reason; their 5p/p is peaked 
well exterior to the DOB, in a zone of inverted density strat- 
ification situated in the outer envelope at logT « 5.2. 

This inversion zone, which is present in both models, 
is responsible for a secondary family of unstable pulsation 
modes uncovered by BOOJUM. These modes, not shown in 
Fig. |3] due to our choice of ordinate range, are found for 
harmonic degrees 1=1 and £ = 2, and are characterized 
by growth rates well in excess of 10 -2 and energy densi- 
ties strong l y con centrated in the inversion zone. Following 
ISaio et aTT il998l) . we use this trapping property to classify 
the modes as non-radial strange modes. We note that the 
iron opacity bump is responsible for both the density in- 
versions and the instability of the strange modes trapped 
within them. Since the strange modes are not related to the 
DOB. we shall not discuss them further. 



5 DISCUSSION 

In the preceding section, we demonstrate that low-degree, 
intermediate radial-order g modes are unstable in a pair of 
stellar models representative of WN stars, due to a re mech- 
anism operating on the deep opacity bump. This finding 
supports the central hypothesis of the paper advanced in 
Sec. Q and likewise confirms WN stars as the high-mass 
cousins of GW Vir stars, exhibiting a relationship parallel 
to that between the high-mass j3 Cephei and SPB stars, and 
the low-mass sdB pulsators. 

However, there are caveats to our results. Beyond the 
various (reasonable) simplifications we adopt in construct- 
ing our stellar models (Sec. and conducting the stability 
analysis (Sec.|HJ, we have neglected the effects of the strong, 
radiatively driven winds that are characteristic of all W-R 
stars. The stellar models incorporate only the evolutionary 
effects of wind-originated mass loss, and the stability analy- 
sis ignore s the radial outfl ow at the outer boundary. As dis- 
cussed bv lCranmeJ dl99rJ) . an outflow can allow ordinarily- 
trapped p ulsation modes to propagate outward through the 
boundary. iTownsendl <l200Cf) demonstrated that this leakage 
of wave energy can result in an appreciable damping of pul- 
sations. If the unstable g modes found herein are to remain 
unstable, the rate of energy input from the re mechanism 
must exceed the damping due to leakage. 

A thorough evaluation of this issue may have to await 
further theoretical advances in understanding how outflows 
influence pulsation. In the meantime, we discuss an obser- 
vational development of particular relevance to the present 
work, one that poses some interesting questions regarding 
the role played by pulsation (whether leaking or not) in 
establishing structure and modulation in W-R winds. Al- 
though searches for periodic light variations in WN8 stars - 
the subtype most characterized by variability - have pr oven 
largely unsuccessful (see, e.g.. iMarchenko et al.lll99a . and 
references therein) , recent M OST observati o ns of the WN8 
star HD 177230 (WR 123) bv lLefevre et alJ d2005l) have de- 
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Figure 3. Differential work functions, plotted as a function of 
envelope temperature, for the 1 = 1 modes of the two stellar 
models having the largest growth rates. The differential is with 
respect to the dimensionless radius x = r/R*, and the work W is 
expressed in units of the total energy of pulsation Eyy ■ 



Table 2. Properties of the unstable £ = 1 and £ = 2 g modes of 
the "WNL" and "WNE" stellar models. The subscripts m i n and 
max indicate the minimum and maximum values of the relevant 
quantities, respectively. 



model 


£ 


-Pmin/hr 


^max/hf ^min ^raax 


?7max 




"WNL" 


1 


10.8 


20.6 -20 -10 


6.67 x 10- 


5 


"WNL" 


2 




No unstable modes 






"WNE" 


1 


3.83 


12.4 -14 -4 


1.67 x 10- 


-3 


"WNE" 


2 


2.81 


8.34 -15 -5 


1.12 x 10~ 


-3 



tected an unambiguous 9.8 hr-periodic signal in the star's 
light curve. 

Pulsation appears a promising candidate for the ori- 
gin of the star's variation. W-R stars have been known 
for some time to b e unstable toward strange modes (e.g., 
iGlatzel et al"lll993Tl : however, there is debate over whether 
such modes are c ompatible with the o bservations of WR 123. 
On the one hand. lLefevre et alJ i2005T) themselves argue that 
the periods typical to strange modes, being on the order of 
the dynamical timescale rd yn , are far too short to match 
the 9.8 hr period detected by MOST. On the other hand, 
iDorfi et all J2006I) reason that as a WN8 star WR 123 is ex- 
pected to have a significantly larger radius, and hence longer 
strange-mode period s, than the helium-star models on which 
iLefevre et alJ i2005|) based their co nclusions. In support of 
their argument. IDorfi et alJ i2006t) present a WN8 model 
having R, = 15.4Rq that exhibits unstable radial strange 
modes with periods on the order of a half day, seemingly 
compatible with the observations. 

Here, there may be a problem. The large radii of the 
models studied bv lDorfi et alJ ( 2006) are, we believe, linked 
to the authors' ab initio assumption of an envelope hydro- 
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gen abundance X = 0.35. Although representative of many 
WN8 stars, this value is inappropriate to WR 123, which is 
well-kno wn to be largely deple ted of hydrogen (X < 0.005; 
see, e.g., ICrowther etaflll995l) . It seems probable that the 
absence of hydrogen in the star is a consequence of it having 
already shed any residual hydrogen-rich envelope via wind 
mass loss. Accompanying the loss of this envelope should be 
a substantial reduction in the star's radius, as can be seen 
by comparing our "WNL" and "WNE" models in Table Q 

Based on this reasoning, hydrogen-poor WR 123 may 
be charact erized by an appre ciably smaller radius than as- 
sumed by iDorfi et alJ J2006), resulti ng once again in th e 
period mismatch noted originally bv iLefevre et all (12005^ . 
This leads us to advance g modes excited by the DOB as an 
alternative explanation for the star's variability. As Fig. [5] 
illustrates, the 9.8 hr period is covered by the unstable 1=1 
g modes of the hydrogen-poor "WNE" model. While we de- 
scribe this model as 'early- type', we once again emphasize 
that the label is to be understood in an evolutionary rather 
than spectroscopic sense. 

The detectable presence of g modes in WR 123 could 
offer unparalleled insights into the star's internal structure; 
these modes penetrate down to the boundary of the convec- 
tive core, and are therefore particularly well suited to aster- 
oseismological analyses. However, it is unclear how best to 
proceed with the necessary initial step of testing our inter- 
pretation of the observations. Pulsation in g modes is often 
diagnosed through time-series monitoring of spectroscopic 
line profile variations (lpv); the predominantl y horizontal 
veloc ity fields generated by these modes (e.g., lUnno et alJ 
1989) tend to concentrate variability in the wings of line 
profiles. Unfortunately, this approach is not feasible in the 
case of W-R stars, because the supersonic wind washes out 
any lpv arising from (subsonic) pulsational velocity fields. It 
remains to be seen whether other approaches can be devised 
that are able to discriminate g modes from other potential 
sources of variability. 

Beyond the specific case of WR 123, the discovery of 
a new class of Wolf-Rayet instability is of great interest in 
itself. Future investigation can now focus on exploring the 
characteristics and extent of the DOB instability in greater 
detail. Specific topics worthy of attention include establish- 
ing red and blue edges in the HR diagram, and determin- 
ing the sensitivity of the instability against changes in in- 
put physics (e.g., elemental diffusion, mass-loss rates, etc.). 
As we have already noted, a proper treatment of pulsation 
at the outflowing outer boundary is also needed. An inte- 
gral part of this treatment must be a better understanding 
of the reverse side of the pulsation-wind interaction: how 
can pulsation at the wind base lead to the modulations ob- 
served in Wolf-Rayet stars? Initial steps toward addre ssing 
this question were made bv lCranmer fc Owockil lll996t) . but 
there remains much work to be done. 
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